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Variational Monte Carlo study of spin polarization stability of fractional quantum 
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We compare ground state energies by variational Monte Carlo of the spin unpolarized Halperin 
331 and the spin polarized Moore-Read (MR) Pfafnan fractional quantum Hall states at half-filling 
of the lowest Landau level (LLL) and the second Landau level (SLL) as a function of small deviations 
around the Coulomb point via the finite thickness effect and direct alterations to the the first two 
Haldane pseudopotentials. In the comparison we find that in the LLL, either the 331 state or the MR 
Pfafnan may be lower in energy depending on the deviations. In the SLL, however, the MR Pfafnan 
is consistently lower in energy except for large deviations. These results suggest that even under 
moderate deviations in the interaction potential (through various physical processes such as finite 
thickness, Landau level mixing, etc.), the MR Pfaffian description is more energetically favorable 
than the Halperin 331 state in the half-filled SLL (i.e. v = 5/2), consistent with recent experimental 
investigations. 
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I. INTRODUCTION 



The direct experimental observation of quasi-particles 
with non-Abelian anyonic statistics is an active area 
of research that may lead to the construction of a 
fault-tolerant topological quantum computei^^!. There 
are several condensed matter systems under study ex- 
pected to exhibit non-Abelian excitations, including one- 
dimensional nano-wires adjoining s-wave superconduc- 
tor^, cold atomic systems with spin-orbit interactional, 
Px + ipy superconductors^, topological insulators on 
superconducting substrates^, thin semiconductor films 
in contact with ordinary superconductors 7 -' and others. 
Perhaps the most discussed system where non-Abelian 
anyons are proposed to exislP*^is in the fractional quan- 
tum Hall effectP (FQHE) in the half-filled second Lan- 
dau leveP^l, i.e., v = 5/2, where strongly interacting elec- 
trons in two-dimensions (2D) are under the influence of 
a strong perpendicular magnetic field. Even though the 
FQHE at v = 5/2 was first observed in 1987^ it remains 
as somewhat of an enigma since the physics of this FQH 
state differs considerably from the well-known Laughlin 13 
and composite fermior l 1 ' 15 ' family of fractional quantum 
Hall states and is the only even denominator fraction 
to be observed in mono-layer quantum Hall systems so 
far. At present, the leading theoretical description for 
v = 5/2, is the spin polarized Moore-Read (MR) Pfaf- 
fian wavefunction 8 introduced in 1991, whose excitations 
have been th eoretica lly shown to have non-Abelian any- 
onic statistic s 1 8 * 16 * 17 : Despite this prediction and consid- 
erable numerical evidence 18 supporting the MR Pfaffian 
description, non-Abelian anyons have yet to be directly 
observed experimentally in quantum Hall systems, al- 
though there are recent tantalizing experimental results 
which are interpreted as being consist ent w ith the 5/2 
FQH quasiparticles being non-Abeliar l 19 * 20 l But there 



exist alternative descriptions for the even-denominator 
FQHE that do not require non-Abelian statistics, and one 
of the leading alternatives is the unpolarized Halperin 331 
stated Part of the Halperin (m,m,n) family of Abelian 
fractional quantum Hall states^, the 331 state is ob- 
served in half-filled quantum hall bi-layers where even 
denominator fractions are commonly observed and well 
understood as belonging to 331 type multicomponent 
Abelian even-denominator FQHE statea^HSlD Resolving 
the true nature of the v = 5/2 (i.e. whether the excita- 
tions are Abelian or non-Abelian) would help determine if 
the fractional quantum Hall effect can potentially serve as 
a testbed for topological quantum computing. For exam- 
ple, if the 331 type state can be convincingly ruled out for 
the 5/2 FQHE, this will provide further indirect evidence 
for the 5/2 FQHE as being the non-Abelian Moore-Read 
type state. 

Determining the spin order (i.e., whether the actual 
5/2 FQH state is spin polarized or not) could possibly 
help resolve the nature of the v = 5/2 state. Since the 
MR Pfaffian (Halperin 331) is, by construction, spin po- 
larized 2 - 7 (unpolarized), determining the spin polariza- 
tion of the state can unequivocally rule-out either the MR 
Pfaffian or the unpolarized 331 description for v = 5/2. 
However the spin-order of v = 5/2 has a long history of 
inconclusive and seemingly contradictory studies. One 
of the earliest experiments on this matter, for example, 
reported the collapse of the v — 5/2 FQHE state in a 
tilted magnetic field, suggesting an unpolarized FQHE 
stated But this experiment was followed by a seminal 
numerical study by Morf 18 which showed that even in the 
limit of zero Zeeman energy, a polarized state has lower 
ground state energy per particle than a fully unpolar- 
ized state at v = 5/2 . Sub sequent theoretical studies re- 
ported similar resulta^^Sl anc ] ft j s now understood that 
the collapse in the original "tilt" experiments was likely 



due to the appearance of compressible "striped" phases^ 
or from orbital coupling of the in-plane magnetic field 32 . 
To avoid these issues, Liu et alP2 performed an "altered" 
tilt experiment where the in-plane magnetic field is kept 
small by varying the charge density and showed that the 
v = 5/2 state remains stable for relatively large tilt an- 
gles. Recent rather impressive experiments performed 
by Tiemann et alP^and Stern et alP" provide strong ev- 
idence supporting a fully spin polarized 5/2 FQHE state. 
But a notable exception is an earlier experiment by Stern 
et alP3 using photoluminescence spectroscopy to probe 
polarization, suggesting that the v = 5/2 state is actu- 
ally unpolarized. A similar result was also seen in a study 
by Rhone et alW* where spin-order is probed via resonant 
light scattering. It is possible that the signatures seen in 
these optical experiments 36 3 ' are due to local spin order 
near a charge impurit y 38 * 39 ^ and are, therefore, inconclu- 
sive, but such results still leave some doubt to the actual 
nature of spin-order in this stated The most convinc- 
ing measurement to date are those by Tiemann et alPS 
and Stern et al.PH which do indicate a spin-polarized 
5/2 FQHE. Our goal in the current work is to provide 
a reasonably complete study of spin-polarization com- 
paring the candidate states MR Pfaffian versus 331 with 
respect to the 5/2 FQHE using direct numerical tech- 
niques in reasonably realistic theoretical models to see if 
one or the other can be ruled out purely on the basis of 
numerical studies. 

An alternative explanation that may resolve the seem- 
ingly conflicting results on spin-order that we explore 
here is the possibility that the v = 5/2 FQHE state is 
not unique^. In other words there may be more than 
one incompressible FQHE state, each with different spin 
polarization, that satisfy conditions to be experimentally 
observed at v = 5/2. Whichever FQHE state is observed 
in an experimental sample may depend on certain de- 
tails of the sample or experimental parameters because 
the competing polarized and unpolarized states may be 
relatively close in energy so that one or the other gets 
stabilized by the minute details of the sample parame- 
ters. 

In this work, we explore this possibility by compar- 
ing the energy of two trial FQHE wavefunctions with 
respect to small deviations in the system Hamiltonian: 
the spin polarized Moore Read (MR) Pfaffian and the 
unpolarized Halperin 331 state. To examine the possibil- 
ity of a phase change between spin polarized and unpo- 
larized incompressible FQHE states within the v = 5/2 
plateau, we focus entirely on the spin polarized MR Pfaf- 
fian and the spin unpolarized Halperin 331 state and ex- 
amine how their respective ground state energies change 
with respect to small alterations to the effective inter- 
action potential. By appealing to the variational the- 
orem, whichever state has the lowest ground state en- 
ergy (strictly speaking, the lowest free energy) provides 
a better physical description-of course, the variational 
theorem cannot rule out the possibility of a lower energy 
ground state that we are not investigating (in fact, this 



is the likely scenario in the nu=l/2 state in the lowest 
Landau level where experimentally an FQHE has never 
been observed indicating that some kind of a compress- 
ible non-FQH state is likely to be lower in energy than 
either the MR Pfaffian or the Halperin 331 state). To 
achieve this goal we alter the 2D Coulomb potential in 
two ways: i) through the finite thickness effect, and ii) 
by directly perturbing the first two Haldane pseudopo- 
tentials (see below) . 

In the finite thickness effect, the non-zero thickness of 
the quasi-2D electron system provides an effective poten- 
tial slightly modified from the purely ideal 2D Coulomb 
potential, which we refer to as the "Coulomb point" in 
this work. Since FQHE states have been shown to be sen- 
sitive to this finite thickness effect and the thickness is 
expected to vary for different experimental samples , the 
finite thickness effect is a natural area to investigat e! 32 * 42 ! 
(in these references the finite-thickness effect was consid- 
ered in the context of spin polarized electrons via exact 
diagonalization) . Although our results do not suggest a 
direct quantum phase transition between the MR Pfaf- 
fian and the Halperin 331 in the second Landau Level 
(SLL) induced by tuning the finite thickness, the energy 
difference between the two states decreases with sample 
thickness, implying that other perturbations may make 
either state energetically favorable for very deep wells; 
some of these other perturbations could be Landau level 
mixing, disorder, effect of nearby gates, self-consistency 
of the confining potential itself, etc., which are not con- 
sidered in our work since they are far too sample-specific 
to be treated theoretically at this stage.. 

Our second approach to altering the 2D Coulomb po- 
tential - directly perturbing the first two Haldane pseu- 
dopotentials33 - provides us a general theoretical probe 
that can identify areas of interest which may be reached 
experimentally via other effects. The Haldane pseudopo- 
tentials, V m , parametrize the effective interaction poten- 
tial in terms of the relative angular momenta m between 
two particles, and thus, perturbing the first two pseu- 
dopotentials (i.e. m = 1,2) alters the short range interac- 
tions between electrons (note that in this study, we leave 
the m = term fixed). The method of altering pseu- 
dopoentials is a common approach taken in FQHE exact 
diagonalization numerical studies aimed at probing sensi- 
tivities to different moments in the interaction strength, 
however, this approach has not yet been attempted in 
VMC studies to the best of our knowledge. We note that 
our two alternative ways of introducing 'small deviations' 
or tuning away from the pure Coulomb point (realistic fi- 
nite thickness effect and varying the lowest pseudopoten- 
tials) are complementary theoretical methods of tuning 
the system Hamiltonian since the finite thickness cor- 
rection modifies all the Coulomb pseudopotentials in a 
complex manner which cannot simply be simulated by 
changing the two lowest pseudopotentials. 

We also examine the lowest Landau level (LLL) (i.e. 
v = 1/2) in addition to the second Landau level (SLL) 
(i.e. v — 5/2) for comparison. In the LLL, no in- 



compressible even-denominator FQHE has been exper- 
imentally observed in mono-layer systems to date, but 
there are several theoretical proposals to engineer cer- 
tain experimental conditions in such a way that even- 
denominator states are energetically favorable^. In these 
cases the Halperin 331 state and the MR Pfaffian are 
likely possibilities. In our study of the LLL, we find that 
either state may be energetically favorable in the LLL 
depending on the pseudopotential deviations. This is in 
contrast to our results in the SLL where we find that 
the MR Pfaffian is generally lower in energy than the 
Halperin 331 state for most deviations examined, sug- 
gesting that the MR Pfaffian description is, indeed, bet- 
ter suited for the half-filled SLL, i.e., v — 5/2. This 
is of course also consistent with the most recent experi- 
mental status of the subject where the SLL 5/2 FQHE 
appears to be spin-polarized. Our work, however, indi- 
cates that the corresponding LLL situation is more del- 
icate, and if an incompressible FQHE is ever observed 
in a half-filled LLL monolayer 2D system, it could ei- 
ther be a MR spin-polarized Pfaffian or a Halperin spin- 
unpolarized 331 state. 

We add a theoretical subtlety here which has some- 
times caused some confusion in the literature. The 
Halperin 331 state in general does not obey the full SU(2) 
symmetry (specifically, the so-called Fock condition nec- 
essary for a spin-independent many-body Hamiltonian 
which must conserve the total spin of the system), and 
cannot therefore be a true eigenstate of the single-layer 
Coulomb Hamiltonian since by definition this Hamilto- 
nian obeys the full SU(2) symmetry because the Coulomb 
interaction is spin-independent -the 331 state was origi- 
nally conceived for the double-layer 2D system where the 
Coulomb interaction does indeed depend on the layer in- 
dex and is in general not SU(2) invariant in the layer 
index. This is, however, not a problem for our VMC 
analysis since we are only interested in comparing ener- 
gies between variational ground states (which do not care 
about the symmetry of the Hamiltonian) and are not try- 
ing to obtain the exact theoretical eigenstate of the sys- 
tem (which would be an impossible task any way, even 
the MR Pfaffian can at best be a good variational state 
for the system and by no means the exact eigenstate). 
In the end, the best theory we can hope for is to obtain 
a variational ground state (MR Pfaffian or 331) which is 
adiabatically connected to the exact ground state of the 
experimental system without any intervening quantum 
phase transition so that the spin-polarization status of 
the variational ground state and the exact ground state 
remains the same. Thus, for our purpose, both the MR 
Pfaffian and the 331 are perfectly (and equally) legiti- 
mate variational choices, and which ever has lower VMC 
energy could be construed as the "correct" ground state 
of the system (at least within the narrow, but very rea- 
sonable, restricted variational choice of only two candi- 
date wavefunctions). 

We also mention that all our work leaves out the triv- 
ial Zeeman energy of the system arising from the ap- 



plied magnetic field creating the Landau levels in the 
first place, which helps the spin-polarized state over the 
spin-unpolarized state. Since the applied field is typically 
rather small for the v = 5/2 FQHE, leaving out Zeeman 
energy (which is trivial to include for any given field) is 
probably a reasonable approximation, but it is helpful to 
remember that even if a spin-unpolarized ground state 
arises from our VMC analysis, the Zeeman energy could 
in principle eventually win over, leading to the experi- 
mental state being spin-polarized. The reverse, however, 
is not true, i.e. if the zero-Zeeman splitting situation 
has a (spontaneously symmetry-broken) spin-polarized 
ground state, it is unlikely that finite Zeeman splitting 
will change the ground state to a spin singlet. 



II. VARIATIONAL MONTE CARLO 

EVALUATION OF THE ENERGIES USING 

EFFECTIVE POTENTIALS 

We use variational Monte Carlo (VMC) methods in the 
same spirit as Refs. PHlandl2T)lto estimate the energy per 
particle of the Halperin 331 state and the MR Pfaffian 
state for altered Coulomb potentials in the lowest and 
second Landau level with up to N — 120 electrons and 
extrapolate to the thermodynamic limit (1/-/V — > 0). To 
examine the finite thickness effect, we use the potential 
derived for the infinite square well potential in a planar 
system which is given, in Fourier space, bjES 
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where I is the magnetic length given by ^Jhc/eB and 
d is the thickness. There are other finite thickness po- 
tentials that we can explore^"^, however we expect the 
infinite square well potential to be sufficient in obtaining 
a qualitative comparison between the two stateiP^Ml 

We also examine the effect of directly perturbing the 
Haldane pseudopotentials^ V m for the Coulomb poten- 
tial in LLL and SLL. In particular, we examine the effect 
of the perturbations Vi -> V\ + AVi and V 2 -» V 2 + AV 2 
for pseudopotentials derived from the Coulomb potential 
(i.e. V(k) = 1/fc). In the planar geometry, the Haldane 
pseudopotentials for the nth Landau level are given in 
terms of the Fourier transform of the effective interac- 
tion potential V Q s(k) by 
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where L n (x) are Laguerre polynomials. Note that we 
are using planar pseudopotentials due to two primary 
reasons: (1) it is simpler to use planar Vm s and (2) 
for such large systems N ~ 100 examined in our studies 
there is very little difference between spherical and planar 
pseudopotentials. 



In order to use VMC methods to estimate wavefunc- 
tion energies, we require an effective potential in real 
space, V e g(r), such that the application of Eq. p| re- 
sults in our perturbed pseudopotentials, V m , on the left- 
hand-side. The immediate difficulty we run into is that 
there is no clear procedure to invert Eq. pi) to obtain 
V c g(r) for arbitrary V^-it is a rather non-unique one- 
to-many mapping. Also, even in the unperturbed case, 
estimating energies in the SLL is not straight-forward 
since most FQHE trial wavefunctions under study do not 
have a closed-form expression in the SLL. To get around 
these difficulties, we chose a variable effective potential 
with fitting parameters, c;, and set these parameters such 
that the result of applying Eq. (TsJ) on the effective po- 
tential very closely matches the perturbed pseudopoten- 
tials. And when examining the SLL, we "simulate" the 
SLL in the LLL by fitting the effective potential within 
the LLL to the perturbed SLL pseudopotentials^, that 
is, we project the SLL into the LLL. Several forms for 
the effective potential hav e been u sed for previous Monte 
Carlo studies of the FQHFpS 49 | 50 [ For our study, we use 
the following form for the effective potential (in units of 
e 2 /d): 
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We choose this form because for large enough M, the po- 
tential fits both even and odd pseudopotentials to a rea- 
sonable degree - only odd pseudopotentials are impor- 
tant when fully polarized or spinless wavefunctions are 
under investigation - and the fits to V m for large m (i.e. 
?7i > M) are generally consistent across different pertur- 
bations, AVi and AV2, allowing us to make fair compar- 
isons between different perturbations. In choosing the 
number of terms, M, in the effective potential, there is a 
trade-off between tighter fits to the pseudopotentials for 
larger M and ease with which the Monte Carlo converges 
- the addition of terms in Eq. pi) leads to an oscillatory 
potential that takes, in general, more iterations to reach 
Monte Carlo convergence. For our study, we use M = 6. 
As an example, we show in Fig. [I] perturbed pseudopo- 
tentials V m for AVi = -0.06 and AV 2 = 0.02 in the SLL 
and the corresponding fitted pseudopotentials resulting 
from a non-linear least squares fit of Eq. (I3j to V m via 
Eq. pi). It is worth noting that the Vm's calculated from 
the effective potential for m > M = 6 are very good 
approximations to the actual values and only differ at 
the level of a fraction of a percentage point (~ 0.6% on 
average) 

Throughout this work we make use of the spherical 
geometry where electrons are confined to a 2D spherical 
surface of radius R with a magnetic mo nopole of mag- 
netic charge Q at the center of the sphereHEU The ra- 
dius of the sphere is determined by the magnetic charge 
Q: R 2 = Q. The magnetic charge for a quantum Hall 
state with N electrons at filling factor v is given by 
2Q = N/v + x where \ 1S the topological shift 52 and 
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FIG. 1. a) Perturbed pseudopotentials V m for AVi = —0.06 
and AVa = 0.02 in the SLL (the Coulomb point is the 
Coulomb interaction in the SLL) and the corresponding fitted 
pseudopotentials from Eq. pi). Note that any difference be- 
tween the pseudopotentials calculated from the real space ef- 
fective potential corresponding to the deviated pseudopoten- 
tials and the deviated pseudopotentials themselves is smaller 
than the symbols on the figure, b) The resulting V c g(r) com- 
pared to the Coulomb potential V(r) — 1/r 



depends on the FQHE state under investigation. Lastly, 
the distance between two electrons on the spherical sur- 
face is taken to be the cord distance. 

In the spherical geometry, the unpolarized Halperin 
331 and polarized MR Pfaffian states are given by 
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sin(0j/2) exp(— z0i/2), and uf = m <8> \a) and 
vf = Vi ® \a), \<t) is the spin ket, A is the anti- 
symmetrization operator, and Pf[M] is the Pfafhan of 
the matrix Mij = (ujv, — vju'A^ 1 . The magnetic 
charge for both states is given by 2Q = 2N — 3. 

To evaluate the energy of some wavefunction -0 via 
variational Monte Carlo we calculate the Hamiltonian ex- 
pectation value 
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where fl = (9, <j)) is the particle position on the sphere and 
H = J2i<j V (\ri — r j I ) . When calculating the energy 

of V'Pfaff we make use of the identity, det M = |Pf [M] \ . 
Before presenting the effects of finite-thickness and di- 
rectly perturbing Haldane pseudopotentials we briefly 
discuss the background energy. It is assumed that there is 
a uniform distribution of positive charge on the spherical 
surface so that the total energy is negative and the elec- 
tron's state represents a stable phase of matter. That is, 
we place N positive charges on the surface of the sphere 
and calculate the interaction energy between an electron 
and the background E e i^bg and the interaction energy of 
the background with itself Eb g -bg- For a pure Coulomb 
interaction this works out to be 
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ber that the radius of the sphere is R -- \[Q. Now, 
strictly speaking, this energy comes about by doing a 
rather trivial integral over the surface of the sphere with 
the distance between particles defined as the cord dis- 
tance instead of the arc distance-in the thermodynamic 
limit both choices are equivalent. 

For our calculations it is a little bit more subtle. We 
are considering electrons projected into the LLL with ef- 
fective potentials that take into account finite thickness, 
electrons completely confined to the SLL, and potentials 
produced by small deviations away from the Coulomb 
point through the direct manipulation of V± and V%. Wc 
find this effective potential through Eq. (p| and to get 
the proper background energy we calculate 
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where r(6) = 2i?sin(#/2). While there is no deep physics 
hidden in the background energy, it is needed to ensure 
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FIG. 2. Energy per particle (in units of e 2 /el) as a function 
of finite thickness, d (in units of the magnetic length, I) for the 
MR Pfaffian (red) and the Halperin 331 state (green) in the 
LLL (a) and the SLL (b). The energy per particle increases 
with thickness and the difference in energy between the two 
states decreases. However, in the LLL, the Halperin 331 is 
always lower in energy than the MR Pfaffian while, in the 
SLL, the opposite is true. In other words, finite thickness 
alone does not apparently drive a spin order transition. 



that the ground state energy per particle has a well- 
defined and finite thermodynamic limit. Further, since 
we are comparing two ground state energies, this back- 
ground energy cancels out in a sense. 



III. GROUND STATE ENERGIES FOR 

EFFECTIVE POTENTIALS IN THE LLL AND 

SLL 

As a thorough numerical check, we first list our results 
for the MR Pfaffian and the Halperin 331 state at the 
Coulomb point with no perturbation in Table [I] and find 
that we are in agreement with Refs. Fill and |2T)1 We note 
that other wavefunctions are also possible, for example, 
a Composite Fermion fermi se a (CFF S), in a polarized 
(P) or unpolarized (UP) varietjl^ 1 ^^ 1 and the Haldane- 
Rezayi singlet stated (HR), all whose energies in the LLL 
and SLL are listed in Table [I] In the LLL, it is clear that 
the lowest energy state is the unpolarized CFFS which is 
a gapless state that does not yield the FQHE, although, 
the energy of the 331 state is very close. However, we 



-0.4 
-042 
-0.44 - 
-0.46 - 
-0.48 - 

-0.5- 

-0.52 - 
0.1 



JMR Pfaffian 
]331 




AV, 



AV, 



0.03 


H ' ' ' -\ 










o nfi 








0.04 


_„-^'--~'"' 






nn? 


^, ** ,- "~ 








^ ^ .— ■"■ 




-s = :=5 " 


0.02 


- ^ ^" 




«•" 


0.04 


- 


0.06 


- 


0.08 


- 



-0.08 -0.06 -0.04 -0.02 0.02 0.04 0.06 0.08 0.1 

AV, 



I 


0.04 
0.03 
0.02=- 
0.01^ 




0.03 
0.06 
0.04 
0.02 


b) 


MR Pfaffian 












0.08 
0.06 
0.04 
0.02 


• I 


> 

<l 







fr 

























-0.01 




-0.02 








■ 






-0.02 


j-0.02 
H-0.03 
H-0.04 
^-0.05 




-0.04 

-0.06 

-0.03 

-0.1 
-0 






R 


^^ffl 






-0.04 
-0.06 
-ona 


1 -0.08 


-0.06 -0.04 -0.02 


0.02 


0.04 0.06 0.08 0.1 












AV, 















FIG. 3. Comparison of Halperin 331 and MR Pfaffian 

in the LLL. a) Energy per particle as a function perturba- 
tion strengths AVi and AV2. b) Energy difference E A = 
-Bpfaff — £331 as function of AVi and AV2. Regions where 
either Halperin 331 (E A > 0) or MR Pfaffian (E A < 0) is 
energetically favorable are denoted. The intermediate region 
denotes area where the energies are within numerical uncer- 
tainty of each other. The star designates the Coulomb point 
for reference. The statistical uncertainty in the energies is 
not indicated on these contour plots for ease of presentation. 
However, it is similar in magnitude to what is presented in 
Fig. [2] but the qualitative effects of the uncertainty is indi- 
cated by the "intermediate" regime where both energies are 
within statistical uncertainty of each other. 



know experimentally 55 that no FQHE has yet been ob- 
served in single layer systems at v = 1/2. In contrast, at 
the Coulomb point in the SLL, the lowest energy state in 
Table | is the MR Pfaffian and, in fact, the MR Pfaffian 
has been routinely experimentally observed at v = 5/2 
albeit at low temperatures and in very high-quality sam- 
ples, indicating that the v = 5/2 FQHE is rather fragile 
with a very small gap with possibly competing states 
with comparable energetics nearby in phase space. 

Since the purpose of this work is to investigate the 
spin polarization of the half- filled lowest and second Lan- 
dau levels via VMC, we will focus exclusively on the 
Halperin 331 versus the MR Pfaffian wavcf unctions. A 
full investigation including all possible ansatz and more 
realistic effective potentials that include finite thickness 
and Landau level mixing 5 ^ is beyond the scope of this 
work and will have to await future works. Our work is in 



FIG. 4. Comparison of Halperin 331 and MR Pfaffian 

in the SLL. a) Energy per particle as a function perturba- 
tion strengths AVi and AV2. b) Energy difference E A = 
Spfaff — £331 as function of AVi and AV2. Regions where 
either Halperin 331 {E A > 0) or MR Pfaffian (E A < 0) is 
energetically favorable are denoted. The intermediate region 
denotes area where the energies are within numerical uncer- 
tainty of each other. The star designates the Coulomb point 
for reference. 



the spirit of a restricted variational study which makes 
sense for this problem since the two candidates we use 
(i.e. MR Pfaffian and Halperin 331) are essentially the 
'only game in town' for incompressible even-denominator 
FQHE states in single-layer systems. In fact, part of the 
motivation of this work is to establish the feasibility of 
this sort of VMC investigation of Hamiltonians described 
by effective potentials. 

Next, we examine the finite thickness effect and how 
it changes the expected ground state energy of the MR 
Pfaffian and the Halperin 331 state. In Fig.[2Fa) we show 
the numerically calculated energies in units of e 2 /el, as a 
function of thickness (in units of the magnetic length, I) 
in the LLL. The Halperin 331 state is consistently lower 
in energy but the gap between the energies of the MR 
Pfaffian and the Halperin 331 state decreases with in- 
creasing thickness. Similar results are seen in Fig. [21(b) 
for the SLL where, by contrast the MR Pfaffian is en- 
ergetically favorable, but the energy difference between 
the two states again decreases with thickness. Part of this 
likely stems from the fact that the overall energy scale is 
shrinking due to the finite thickness effect. We emphasize 



TABLE I. Energy calculated via VMC of various wavefunc- 
tions. Our results for -E331 and Ep^a agree with those of Di- 
mov et alW* at the Coulomb point in the LLL and SLL. The 
results listed below for the polarized and unpolarized CFFS 
and the Haldane-Rezayi singlet state are given by Park et 
ffl £|44] rp^g ] owes £ ener gy state at the Coulomb point in the 
LLL and SLL is the unpolarized CFFS and MR Pfaffian, re- 
spectively (both are indicated in bold). 



Energy [e 2 /{d)] 



LLL 



SLL 



£331 

-Epfaff 

-E'CFFS(P) 

-ECFFS(UP) 

-Ehr 



-0.4631(3) 

-0.4573(3) 

-0.46557(6^ 

-0.46953(7)P 

-0.3147(3^ 



-0.329(3) 
-0.361(2) 

-0.3492(5^ 
-0.2952(3^ 
-0.303(3JS1 



that the deceptive simplicity of Fig. [2] hides the complex- 
ity of the calculation. For each point on both (rather the 
four) curves the following procedure was carried out: (1) 
for each value of thickness d the pseudopotentials were 
calculated, (2) a real space effective potential V c s(r) was 
found from these pseudopotentials, (2) many VMC eval- 
uations of the energy of either \I>3 31 or Epfaff for N elec- 
trons were carried out, and finally (4) these energies were 
extrapolated to the thermodynamic limit (1/N — > 0) to 
generate a single point. The numerical work in producing 
Fig. [2] is very exhaustive. 

We now examine the effect of directly perturbing the 
first two Haldane pseudopotentials in the LLL and the 
SLL. Fig. pJFa) gives the energy estimates for the MR 
Pfaffian and the Halperin 331 state with respect to AVi 
and AV2 in the LLL. Here we see that the energy of the 
MR Pfaffian is more sensitive to AV m compared to the 
Halperin 331 state. Also the MR Pfaffian decreases in en- 
ergy with decreasing AVi or increasing AV2, whereas the 
Halperin 331 energy shows opposing trends. Figure^ (b) 
shows the energy difference between the MR Pfaffian and 
the Halperin 331 state and an estimated phase diagram- 
the phase is determined by whichever wavefunction has 
the lowest ground state energy per particle. The "inter- 
mediate" phase indicates where the energy differences of 
the two states are within numerical uncertainty of the 
VMC. Here we see that the MR Pfaffian can be energeti- 
cally favorable in the LLL for relatively small deviations 
from the Coulomb point for AV2 > although the 331 
state is lower in energy at the LLL Coulomb point. 

Results for perturbations about the SLL Coulomb 
point are given in Fig. Hi a). Here we see the energy ver- 
sus AVi and AV2 has similar trends as was found in the 
LLL , but the MR Pfaffian is consistently lower in energy 
in the SLL for the majority of perturbations under inves- 
tigation. Fig. Bb) gives the difference in energy and an 
estimated phase boundary between the MR Pfaffian and 
the Halperin 331. Again, the "intermediate" area indi- 
cates where the energies are within numerical uncertainty 
of each other. Unlike in the LLL case, the MR Pfaffian 



is generally favored for any perturbation in the SLL. In 
the region where the Halperin 331 state is favorable in 
the SLL, the perturbations result in V 2 < V\, which is a 
qualitative feature of the Coulomb point in the LLL. If 
the effect of V m> 2 are minimal, then we can argue that 
this region is qualitatively similar to the Coulomb point 
of the LLL and therefore, the Halperin 331 state is ener- 
getically favorable in this region given the results on the 
LLL. 



IV. CONCLUSIONS 

Our results show that through the finite thickness ef- 
fect, the energies of the Halperin 331 and the MR Pfaf- 
fian state increase with increasing sample thickness in 
either LL. The energy difference between the two states 
decreases with increasing sample thickness, but there is 
no crossing between the two states in either LL, i.e., finite 
thickness apparently does not drive a spin polarization 
quantum phase transition at least for the situation with 
vanishing Zeeman energy considered in our work. It is 
possible, in fact quite likely for the LLL, that a finite Zee- 
man splitting will induce a transition from the Halperin 
331 to the MR Pfaffian state, but neither may be the true 
ground state in the LLL since no v = 1/2 FQHE has ever 
been observed experimentally. Additionally, our results 
show that the energy of the MR Pfaffian is more sensi- 
tive to changes in the pseudopotentials than the Halperin 
331 state, where the MR Pfaffian energy decreases with 
increasing AV 2 and decreasing AVi. The energy of the 
Halperin 331 state, in contrast, is generally insensitive (in 
comparison to that of the MR Pfaffian) with slight de- 
creases for AVi > and AV 2 < 0, cf. Figs.^a) andga). 
In the LLL, the MR Pfaffian becomes energetically favor- 
able for small increases in V\ above the Coulomb point. 
In the SLL, the Halperin 331 state becomes favorable 
for relatively large deviations from the Coulomb point, 
(i.e. AVi ^ and AV 2 < —0.06), again in the absence 
of the Zeeman energy - given the small difference be- 
tween the two VMC energies, it is quite likely that the 
MR Pfaffian state has lower energy than the Halperin 
331 state for all thickness values and all deviations from 
the Coulomb point once the Zeeman energy is taken into 
account since the Zeeman energy would always prefer the 
spin-polarized state. 

Our study adds to the growing body of 
evidence 1 18 * 30 * 33 " ^ supporting the spin polarized 
MR Pfaffian description for FQHE at v = 5/2. Of 
course, conclusive verification of the MR Pfaffian de- 
scription requires the direct experimental observation 
of non-Abelian anyons. But given the dif ficulty i n con- 
clusively detecting non-Abelian signature d 19 * 20 * 57 !, novel 
experimental techniques will be needed for definitive 
verification. 

Lastly, we emphasize that our study additionally serves 
as a "proof of principle" for VMC studies of various FQH 
systems that are described by effective potentials. Effec- 



tive potentials are needed when considering certain real- 
istic effects such as finite thickness, Landau level mixing, 
higher Landau level FQHE, disorder, etc., or by simply 
artificially manipulating various Haldane pseudopoten- 
tials. Our work establishes that the VMC technique is 
a viable alternative to the exact diagonalization method 
in theoretically studying the ground state properties of 
FQHE including various realistic effects which are often 
hard to implement using the finite size diagonalization 
method. 
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